Lorentz shear modulus of a two-dimensional electron gas at high magnetic field 
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We show that the Lorentz shear modulus - one of the three elastic moduli of a homogeneous 
electron gas in a magnetic field - can be calculated exactly in the limit of high magnetic field (i.e. 
in the lowest Landau level). Its value is ±ftn/4, where n is the two-dimensional electron density 
and the sign is determined by the orientation of the magnetic field. We use this result to refine our 
previous calculations of the dispersion of the collective modes of fractional quantum Hall liquids. 



The collective dynamics of a two-dimensional electron 
gas at high magnetic fields continues to be a subject of 
intense interest and a source of surprises. Witness be 
the recent observation, by inelastic light scattering, of 
a new collective mode in the fractional quantum Hall 
liquid at filling factor v — X/Z} Theoretically, the inter- 
est arises from the fact that the electrons in the lowest 
Landau level (LLL) exhibit an unusual kind of collec- 
tive behavior (for the most recent overview of the theory 
of fractional quantum Hall effect see, e. g., Ref. i). In 
an ordinary two-dimensional quantum liquid collective 
modes arise either as hydrodynamic modes (sustained by 
frequent collisions in local quasi-equilibrium) or as coUi- 
sionless modes sustained by strong self-consistent fields 
(in which case a generalized hydrodynamic description 
is possible). In both cases the frequency of the modes 
tend to zero at long wavelength, reflecting the under- 
lying translational invariance of the system. By con- 
trast, in the incompressible fractional quantum Hall liq- 
uid the collective modes have a finite frequency in the 
long wavelength limit. One mode at the classical cy- 
clotron frequency = ^ (— e and m being the electron 
charge and mass ~ we use SI units throughout) is com- 
pletely expected as a consequence of Kohn's theoren*^. 
The surprise comes from the low-frequency mode, which 
in a fractional quantum Hall liquid (e.g. the v — X/Z 

state) tends to a finite frequency A which scales as 

{£ — \J'-^ being the magnetic length). The gap has long 
been understood as a manifestation of the incompress- 
ibility of the electron liquid in the LLLM, which is also 
responsible for the fractional quantum Hall effect^. (The 
gap closes in compressible states, e.g. the v = 1/2 state.) 

The presence of the gap raises the question whether 
a continuum mechanics formulation (hydrodynam- 
ics/elasticity) for this type of modes is possible at all. 
At first sight such a formulation can only produces gap- 
less modes. Yet, the absence of low-energy excitations or 
other single particle excitations overlapping the energy 
of the collective mode strongly suggests that a collective 
description of the dynamics should be possible. Early at- 
tempts to formulate such a description relied on the ad- 
hoc introduction of a diverging bulk modulus^. Recently, 



it has been found that the problem has a natural formula- 
tion as an elasticity theory in which the elastic constants 
are local in space, but strongly retarded in time^i^ii^. Let 
us briefly review the essential points of this formulation. 

First, the equation of motion for a Fourier component 
of the particle current density at frequency w, in the lin- 
ear approximation and in the absence of external fields 
(other than a magnetic field B in the negative z direc- 
tion) is 

— «mci;j(r, w) — ei?j(r, x z + V- P (r, w) = (1) 

where — e and m are the electron charge and mass respec- 
tively, P (r, ljj) is the Fourier component of the stress ten- 
sor, and V- denotes its divergence. The stress tensor is 
defined as the expectation value of the stress tensor oper- 
ator P, which can be directly derived from the Heisenberg 
equation of motion for the current operator— li^. 

Second, the stress tensor is expressed (in the long- 
wavelength limit) as a linear function of the current den- 
sity 



Py(r,cj) = -Qijki{uj)uki{v,ui) 
(we sum over repeated indices) where 



(2) 



Uki (r, w) = - — [dkji {y,uj)+ dijk (r , w)] , (3) 
znu! 



is the Fourier component of the strain tensor (n being the 
homogeneous density of the electron gas), and Qijud^) 
is the Fourier component of the homogeneous tensor of 
elasticity - a fourth-rank tensor, which is symmetric with 
respect to any of the interchanges {ij) {kl), i <-> j, or 
k ^ I. In two dimensions Qijki{i^) has at most 6 indepen- 
dent components. However, in a two-dimensional homo- 
geneous electron liquid in the presence of a perpendicular 
magnetic field the six components of Qijki (i^) can be ex- 
pressed in terms of only three independent elastic mod- 
uli in the following manner: Qxx.xxii-^) — Qyy,yy{^) — 

K{uj) + Hiuj), Qxy,xy[<^) = K^) , Qxx,yy{i^) = K (uj) - 
fl{Lu), and Qxx,xyiuj) = -Qyy,yx{Lu) = iujA{Lu). A com- 
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pact representation is 

■ ^{'^) / r r r r N 

+ ^ [SikOjl + EjkOil + EilOjk + ejlOik). 

(4) 

where K{ijj), ^{lo) and A(cj) arc the frequency-dependent 
bulk modulus, the shear modulus, and the Lorentz shear 
modulus respectively^'^. The first two moduli, K and 
/i, are familiar from conventional elasticity theoryi^, and 
cannot be calculated without detailed microscopic input 
(notice however that the zero frequency limit of the shear 
modulus must be zero in a liquid state). The Lorentz 
shear modulus A, introduced in Refs. [slfolliol. is a novel 
feature of the system in the presence of a magnetic field. 
Physically it is responsible for the force that squeezes 
together two oppositely directed streamlines. Formally 
it arises as follows: First, rotational invariance about the 
z axis mandates that the tensor of elasticity be invariant 
under the transformation x y, y ^ —x: this enforces 
the identity Qxx,xy{'^) = —Qyy,yx(y) at all frequencies. 
Second, at zero frequency one has an additional reflection 
symmetry which allows us to interchange x and y without 
the minus sign: this implies Qxx,xy{^) = ~Qyy,yxiO) — 0. 
In the absence of a magnetic field this symmetry persists 
even at finite frequency, but the magnetic field breaks it: 
hence A(a;) has a finite non-zero value. 

The main purpose of this Communication is to show 
that the zero frequency limit of the Lorentz shear mod- 
ulus Aq = — lim^j^o A(w) can be calculated exactly and 
has the universal value ±?in/4 (-1- if the magnetic field 
is along the negative z axis, — if it is along the posi- 
tive z axis) when all the electrons are in the lowest Lan- 
dau level. The reason is that Aq can be expressed as a 
Berry curvature of the ground-state wave function with 
respect to a homogeneous change of the metrics, just as 
the Hall conductivity can be expressed as a Berry curva- 
ture with respect to a homogeneous vector potentia l "'^^i^^ . 
This analogy was first noted in a paper by Avron et al^, 
entitled "The viscosity of the quantum Hall liquid" , in 
which Aq was calculated for a non interacting electron 
gas at integral filling factor, and improperly called "vis- 
cosity". Because of the non-dissipative character of the 
dynamics we will continue to refer to Aq as the "Lorentz 
shear modulus". Unlike Avron et ai, who considered 
only noninteracting systems at integral filling factors, we 
will calculate the Lorentz shear modulus for an interact- 
ing electron gas in a partially occupied lowest Landau 
level. We will show that the calculation can be done ex- 
actly due to a special analyticity property which every 
wave function in the LLL enjoys. Armed with the exact 
Lorentz shear modulus we will then return to the equa- 
tion of motion ([1]) and show how the dispersion of the 
collective modes obtained in RefJ^ by treating Aq as a 
fit parameter, gets modified when the exact value of Aq 
is used instead. 

A very useful representation of the stress tensor can be 
given in terms of the derivative of the Hamiltonian with 



respect to a metric tensor 5y (r) [ds^ = Qijdx^dx^) which 
defines a non-euclidean geometry in the plane {gij = 5ij 
is the Euclidean metrics) (see, e. g.^^^). This is the ana- 
logue of defining the current operator as the derivative of 
the Hamiltonian with respect to the vector potentiali^. 
To this end, let us introduce the Hamiltonian of the spin- 
less two-dimensional electron gas in a metric field gij (r) : 

where the sums run over particle indices (n and n'), 

r„ ^ -^g-i (rn)DnW g{rn)9"^ {rn)Dn,jg^Hrn) (6) 
Zm 

is the kinetic energy of the n-th electron, is the inverse 
ofg,j, 

Dns = dr„i + ieAi{rn) (7) 

is the kinetic momentum operator, A(r) is the vector po- 
tential, and d[r, r'] is the length of the geodesies connect- 
ing r to r' in the non-Euclidean plane. We will work in 
the Landau gauge A(r) — {By, 0) with periodic bound- 
ary conditions in the x direction and open boundary con- 
ditions in the y direction. The linear size of the system 
in either direction is L, which we take as our unit of 
length, L ~ 1. We have also set ft = 1. In these units the 
operator D is conveniently expressed as 

Dx = dx + 2TTiNLy, Dy = dy, (8) 

where is the number of magnetic fiux quanta (h/e) 
enclosed by the system. 

The stress tensor is defined as follows 

where the average is taken in the time-dependent quan- 
tum state and at Euclidean metric. With this definition 
of the stress tensor, the equation of motion ^ is ex- 
act. The basic task of the theory is to express the stress 
tensor as a (linear) functional of the current density so 
that Eq. ([1]) can be closed. This is a formidable task, of 
course, but we know that it can be carried out exactly in 
principle^-. 

As a first step in this direction we notice a small de- 
viation from the Euclidean metrics, arising from the in- 
finitesimal deformation r — > r-l-u(r) (u being the integral 
in time of the velocity field j/n), can be represented in 
the form 5*-' = Sij — 2it,y where is the strain tensor. 
Then, it is not difficult to show that Qijki (i^) has the 
microscopic representation^^ 

Q^Jkl{uJ) = Q°;^,^ + m,;Pkl))u, (10) 

where the first term (purely real and independent of fre- 
quency) is the instantaneous derivative of Py with re- 
spect to g'^' and the second term is the Kubo formula for 
the stress-stress response function. 
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It is evident from Eq. (|4|) that 



An = — lim 



Qxxxy (^) 



UJ 



(11) 



Then, using Eq. (jlO[) in combination with the geometric 
definition of the stress tensor Q we arrive, after weU- 
known manipulations analogous to the ones that lead to 
the expression for the Hall conductivit y^^'^^ , to the key 
result 



Aq = 8 X 3m 



QgXX 



dg^y 



(12) 



where |V'o) is the ground-state of the Hamiltonian ([5]) 
in the presence of a homogeneous metrics g''^ . This is 
essentially the formula obtained in Ref. [l3 for what they 
call "viscosity". 

In order to proceed, it is very convenient to 
parametrize the metric tensor as follows 



J 

ffy = — 

T2 



1 



1 



, and g*-' = — — 



(13) 

where t — ti + 1x2 is a complex number defining the 
length and orientation of the y axis of the distorted plane, 
and J is the Jacobian of the coordinate transformation 
that is induced by deformation. The Euclidean metrics 
is recovered by setting J = 1, n = and T2 — 1- The 
derivatives with respect to g^^ evaluated at the Euclidean 
metrics are given by 



d 



QgXX 

d 



dg^y dgy 



1 /_9 d_ 

2 \dT2 ^ dJ 
d _ 1 d 



d 



dgy 



1 f _d_ 

2 Ut2 



d_ 
dJ 



(14) 



Notice that the derivative with respect to g^^ is taken at 
constant g^^ and viceversEfc^^. 

Inserting Eq. ([T^ into Eq. (|12p. and noting that the 
wave function does not depend on J (J enters the homo- 
geneous Hamiltonian as a global scale factor), we get the 
following representation for Aq 



OTl 



dT2/' 



(15) 



T2 = 1 



where the derivatives are calculated at ri 0, 
and J = 1. 

The calculation of Aq is dramatically simplified if we 
assume that the wave function lies within the lowest Lan- 
dau level. First of all, let us define "Landau levels" . It is 
easy to see that the eigenfunctions of the kinetic energy 
operator ^ with a homogeneous metric and the stated 
boundary conditions have the form (up to a normaliza- 
tion constant) 



i^ik{x,y) = e^-'^-e^^^^-S Hi{./2^^^y) (16) 

where y = y + k/Ni^, k is an integer, and Hi is the 1-th 
Hermite polynomial. (The eigenvalue is (/ -I- l/2)u!c/J)- 
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FIG. 1: (Color online) Dispersion of the collective modes at 
v = 1/3 obtained from the solution of Eq. ([1} for (i) the value 
of Ao — 1.4/i/(4n) used in Ref. ilO| (dashed lines) and (ii) the 
exact value from Eq. II 1911 fsolid lines). All other parameters 
are the same as in Ref. [lO. Here Q = q£ and uj is in units of 



The (degenerate) states of the lowest Landau level are 
the ones with I = and k ranging from to — 1. 
From this, we immediately see that any A'^-electron wave 
function that lies entirely in the lowest Landau level must 
have the form 



TV 



■0(ri, rjv) = A{Ti,T2)f{£,i, ...,£,n) ^ 



(17) 



where / is an analytic function of the variables = a;^ -|- 
Tj/i, and therefore also of r = ri -I- 1x2. A(ti,T2) is the 
normalization constant. Unlike /, A depends separately 
on Tl and r2, i.e. it is a non-analytic function of t. [We 
note that the form of Eq. ([T7]) is, in fact, independent of 
boundary conditions.] Hence in general the r-dependence 
of any LLL many-body wave function can be represented 
as '0(ti, T2) = A(ti, r2)$(r), where ^{t) is a holomorphic 
function. Making use of this fact it is possible to show 
that Eq. (fT5|) simplifies as foUows^^ 



h 

I2 



drl 



lnA(Tl,T2) . 



(18) 



(Here we have reinstated the physical units, multiplying 
the result of the dimensionless calculation by the factor 
h/L"^, which had been previously set to 1). 

The problem is finally reduced to calculating deriva- 
tives of the normalization constant. A second major sim- 
plification occurs in the thermodynamic limit, i.e. when 
the number N of electrons tends to infinity in such a way 
that N/Nl — V remains constant. Because tends to 
infinity it is evident that only values of y ~ '^/Nl ~ 
contribute to the normalization integral, and therefore 
the dependence of the analytic factor / on r (through 
^ = a; -|- ry) effectively disappears. Then we immediately 
see that the r-dependence of the normalization constant 

is of the form A{ti,T2) oc t^^^, and use of Eq. yields 



An 



hn 

T 



(19) 
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where n = N/L"^ is the areal density of electrons. We em- 
phasize that this result does not depend on the detailed 
form of the wave function. 

Let us now apply the above universal result to the 
dispersion of collective modes in incompressible liquid 
states. In Ref. [l3 we developed a theory of the elas- 
tic constants in the incompressible quantum Hall liquid, 
based on the premise that the the stress-stress response 
function is dominated by a single pole at the frequency A 
of the q = gap. This assumption uniquely fixes the fre- 
quency dependence of the ordinary shear modulus ^i{uj) 
and the Lorcntz shear modulus 



AoA2 



while the bulk modulus K remains frequency indepen- 
dent for any LLL state. The excitation gap A, the bulk 
modulus K together with /ioo and Aq became the phe- 
nomenological parameters of the theory and A how- 
ever could be determined self-consistently from sum rule 
arguments provided the microscopic ground state wave 
fimction is given^°). 

Inserting the above fJ.{uj), A(cj), and K into Eq. (j4|) 
and this into Eqs. ^ and H]), and solving the result- 
ing dispersion equation we get two collective modes with 
frequencies u!±{q) 



^ = 1 + Moo(/ioo + K)^ - Aog2 



(21) 



± 



Moo(Moo + K) 



AoQ2 



+ (/2^ - A2)Q4 



the dimensionless wave vector and elastic moduli. 

In the special case of a Laughlin wave function with 
filling factor ly = l/M, where M is an odd integer larger 
than 1, we have shown in Ref. [13 that ^ — ^HilMiLLi, 
Combining this with our new result of Eq. (jl9p we get 

AoA _ 



where Q ^ qi, fl^ ^ K ^ and Ao - are 



Hence the only phenomenological param- 
eter that is left at the current stage of the theory is the 
bulk modulus K. Importantly, K is irrelevant for small 
q. Therefore at long wavelength our dispersion relations 
pTjl for the Laughlin states are completely parameter- 
free. 

Fig. [T] shows the dispersion of the collective modes cal- 
culated at = 1/3 with the exact value of Aq in com- 

(20) 

parison with the somewhat larger value of 1.4-^ used in 
Ref. [13 (the bulk modulus remains the same) . It is clear 
that the overall behavior of the curve is still quite satis- 
factory. The reduced values of Aq flattens the dispersion 
of the collective modes at large q giving a shallower roton 
minimum. The changes in the dispersion at small q are 
quite small. 

In summary, we have presented an exact calculation 
of the low-frequency elasticity modulus Aq of interacting 
electrons in the lowest Landau level. The fact that the 
calculation can be done in universal form is a consequence 
of two facts: (i) The existence of a geometric represen- 
tation of Aq, and (ii) The analyticity of lowest Landau 
level wave functions with respect to the complex defor- 
mation parameter r, which describes a variation of the 
metrics. We believe that this general result will be useful 
to further elucidate the dynamics of collective states (not 
necessarily uniform) in the lowest Landau level. 
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